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Abstract 

We report on a search for flavor-changing neutral-currents (FCNC) in the 
production of heavy quarks in deep inelastic Vfj,N and u^N scattering by 
the NuTeV experiment at the Fermilab Tevatron. This measurement, made 
possible by the high-purity NuTeV sign-selected beams, probes for FCNC in 
heavy flavors at the quark level and is uniquely sensitive to neutrino couplings 
of potential FCNC mediators. All searches are consistent with zero, and limits 
on the effective mixing strengths l^^d^j l^dftl^) and \Vsh^ are obtained. 
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I. INTRODUCTION 



Flavor-changing neutral current (FCNC) interactions of c- and 6-quarks appear in a num- 
ber of extensions to the Standard Model (SM) of particle physics including extra quark gen- 
erations [?,?], technicolor [?,?,?,?,?], multiple Higgs sectors (as in supersymmetry) [?,?,?,?], 
left- right symmetric models [?], and leptoquarks [?,?]. Evidence for FCNC effects in the 
heavy quark sector beyond higher order SM processes has not yet been observed. Present 
limits on FCNC result from searches for rare decays of charm [?] and beauty mesons 
[?,?,?], in particular decays of the type D — > t^l''X and B — > where ^ = e or 

H, D = (^D~^ , , Dg^ , B = {B~^,B^,Bg); and X is nothing, a pseudoscalar, or a vector 
meson. The c ^ u transitions are particularly sensitive to new physics since loop level SM- 
FCNC decays are severely suppressed by the Cabibbo-Kobayashi-Maskawa (CKM) matrix. 
While experimental signatures for FCNC in D and B decays are clear, their interpretation 
is ambiguous. Meson decay rates depend on one or more incalculable hadronic form factors. 
In addition, experimentally attractive final states such as D° — > e'^e~ and S° — > ^'^ ^~axe 
helicity-suppressed, which obscures dynamical roles played by particular FCNC models. 

This article presents an alternative search for FCNC processes in the DIS data of the 
NuTeV experiment; where either neutrinos or anti-neutrinos interact with a massive iron 
target. Flavor changing effects can be sought in the reactions 

c ^ fi+X', (1) 
b ^ 1,+X', (2) 
b cX', c i^+X", (3) 

and their charge-conjugates. The experimental signature in the detector is a muon of oppo- 
site lepton number from the beam neutrino. It is possible to isolate this final state because 
NuTeV ran with a high purity sign-selected beam in which the i^n/yn event ratio in neutrino 
mode and the u^/u^ ratio in anti-neutrino mode were 0.8 x 10~^ and 4.8 x 10~^, respectively. 
Because of the semi- inclusive character of the measurement, FCNC effects in neutrino scat- 
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v^N u^cX, 
v^N i^ijbX, 



tering can be probed at the quark, rather than the hadron level. Furthermore, neutrino 
scattering is particularly sensitive to any FCNC process mediated by an intermediate neu- 
tral object that couples more strongly to neutrinos than to charged leptons (e.g., a Z°-like 
coupling) . 

II. EXPERIMENTAL APPARATUS AND BEAM 

The NuTeV (Fermilab-E815) neutrino experiment collected data during the 1996-97 fixed 
target run with the refurbished Lab E neutrino detector and a newly installed Sign-Selected 
Quadrupole Train (SSQT) neutrino beamline. The sign-selection optics of the SSQT pick the 
charge of secondary pions and kaons which determines whether i/^ or i/^ are predominantly 
produced. During NuTeV's run the primary production target received 1.13 x 10^^ and 
1.41 X 10^^ protons-on-target in neutrino and anti-neutrino modes, respectively. The SSQT 
and its performance are described in detail elsewhere [?]. 

The Lab E detector [?], consists of two major parts; a target calorimeter and an iron 
toroid spectrometer. The target calorimeter contains 690 tons of steel sampled at 10 cm inter- 
vals by 84 3mx3m scintillator counters and at 20 cm intervals by 42 3mx3m drift chambers. 
The toroid spectrometer consists of four stations of drift chambers separated by iron toroid 
magnets. Precision hadron and muon calibration beams monitored the calorimeter and 
spectrometer performance throughout the course of data taking. The calorimeter achieves 
a samphng-dominated hadronic energy resolution of oehadI ^had — 2.4% ® S7%/y/ Eh ad 
and an absolute scale uncertainty of 5Ehad/ Ehad — 0.5%. The spectrometer's multiple 
Coulomb scattering dominated muon energy resolution is Ge^/ = 11% and the muon mo- 
mentum scale is known to dE^/ E^j^ — 1.0%. With the selection criteria used in this analysis, 
the muon charge mis-identification probabihty in the spectrometer is 2 x 10~^. This latter 
rate is confirmed by measurement the muon calibration beam. 

III. ANALYSIS PROCEDURE 
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A. Introduction and Data Selection 



The analysis technique consists of comparing the distributions of yvis — 
Eh ad/ {Eh ad + E^) measured in the i/^ and wrong sign muon (WSM) data samples 
to a Monte Carlo(MC) simulation containing all known conventional WSM sources and a 
possible FCNC signal. The FCNC signal peaks at high yyis because the decay muon from 
the heavy flavor hadron is usually much less energetic than the hadron energy produced in 
the NC interaction. The largest backgrounds, from beam impurities, are concentrated at 
low Uvis in ^ij. a-nd distributed evenly across yvis in ^fj. mode due to the respective (1 — y)^ 
and uniform-in-?/ characteristics of the CC interactions of wrong-flavor beam backgrounds. 

Events in the WSM sample must satisfy a number of selection criteria("cuts"). The 
fiducial volume cut requires that event vertices be reconstructed at least 25 cm-Fe (cm of 
iron-equivalent) from the outer edges of the detector in the transverse directions, at least 35 
cm-Fe downstream of the upstream face of the detector, and at least 200 cm-Fe upstream 
of the toroid. Events must possess a hadronic energy of at least 10 GeV, and exactly one 
track (the muon) must be found. The muon is required to be well-reconstructed and to pass 
within the understood regions of the toroid's magnetic field. The muon's energy must be 
between 10 and 150 GeV, and its charge must be consistent with having the opposite lepton 
number as the primary beam component. Requiring that the muon energy reconstructed in 
different longitudinal sections of the toroid agree within 25% of the value measured using the 
full toroid reduces charge mis-identification backgrounds to the 2 x 10~^ level. Finally, for 
the purposes of the final FCNC fit, the reconstructed yvis is required to be larger than 0.5. 
With these cuts there are 207 v-mode and 127 P-mode WSM events remaining in NuTeV's 
nearly 2 million single muon sample. 
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B. Source and Background Simulations 



Conventional WSM sources arise from beam impurities, right-flavor charged current (CC) 
events where the charge of the muon is mis-reconstructed, CC and NC events where a tt or 
K decays in the hadron shower, charged current (CC) charm production where the primary 
muon is not reconstructed or the charm quark is produced via a interaction, and neutral 
current (NC) cc pair production. Single charm CC production and NC cc pair production 
background sources produce broad peaks at high yvis and must be handled with care. 
Table ?? gives the fractional contribution of each background component. The relatively 
large beam impurity background consists of contributions from hadrons (including charm) 
that decay before the sign-selecting dipoles in the SSQT, neutral kaon decays, muon decays, 
decay of hadrons produced by secondary interactions in the SSQT ( "scraping" ) , and from 
decay of wrong-sign pions produced in kaon decays. Table ?? summarizes the relative 
contributions of each beam source. 

A complete GEANT [?] simulation of the SSQT is used to model beam impurities. 
This simulation uses Malensek's [?,?] parameterization for hadron production from the pri- 
mary target. Scraping contributions are modeled by GHEISHA [?]. Production of is 
handled by extending Malensek's charged kaon parameterizations using the quark count- 
ing relation = {?,K~ -\- K^) / A. Charm production is parametrized using available 
data from 800 GeV proton beams [?,?]. GEANT properly handles cascade decays such as 
7r=^7r=^7r'F,7r'F iJ^V^{i'^) and //=^P^(i/^), //=^ e=^P^(i/^)i/e(i^e)- The NuTeV 

detector is hkewise modeled with a GEANT-based hit-level MC simulation. Wrong-sign 
muons generated from the flux simulation are propagated through the detector MC and 
then reconstructed using the same package that is used for data reconstruction. A fast 
parametric MC is also used to compare the high statistics right-sign flux simulation to data 
in i/ju and 1/^. These comparisons showed that the SSQT dipoles required a downward shift 
of -2.5% from their nominal values. The right-sign comparisons after these shifts, are shown 
in Fig. ?? and indicate agreement between predicted flux and data at roughly the 2% level. 
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The high density target-calorimeter suppresses WSM contributions from tt/X decay 
in the hadron shower; their contribution is estimated from a previous measurement of 
/^-production in hadron showers using the same detector [?]. The small charge mis- 
identification contribution is estimated by passing a large sample of simulated events through 
the full detector MC and event reconstruction. 

After impurities, the next largest WSM source comes from CC production of charm in 
which the charm quark decays semi-muonically (dimuon) and its decay muon is picked up 
in the spectrometer while the primary lepton is either an electron or a muon which exits 
from or ranges out in the calorimeter. The v^. beam fraction is 1.9(1.3)% in v (i/)-mode, 
and 22% of the CC charm events which pass WSM cuts originate from a v^. The CC charm 
background is simulated using a leading-order QCD charm production model with produc- 
tion, fragmentation, and charm decay parameters tuned on neutrino dimuon data collected 
by NuTeV [?] and a previous experiment using the same detector [?] . Overall normalization 
of the source is obtained from the measured charm-to-total CC cross section ratio and the 
single muon right-sign data sample. Simulated dimuon events are passed through the full 
GEANT simulation of the detector. Fig. ?? provides a check of the modeling of this source 
through a comparison of the distribution of y'yjs — Eh ad/ {Eh ad + -£'^2), where E^2 is the 
energy of the WSM in the event, between data and MC for dimuon events in which both 
muons are reconstructed by the spectrometer. This distribution should closely mimic the 
expected background to the yvis distribution in the WSM sample. A comparison test 
between data and model yields a value of 19 for 17 degrees of freedom. 

Finally, NC cc production produces a WSM when the c (c) decays semi-muonically in 

(p/j) mode. An excess over other sources at high yvis indicates that this source is present 
in the data; its analysis [?] will appear in a forthcoming publication. For the FCNC search, 
NC charm production is simulated at production level by a gluon fusion model [?] with 
charm mass parameter rUc — 1.70 ± 0.19 GeV/c^ taken from a next-to-leading (NLO) order 
QCD analysis of CC charm production [?] and using the GRV94HO [?] gluon parton 
distribution function (PDF). The NLO charm mass is used because it is influenced in part 
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by contributions from VT— gluon fusion diagrams similar to the Z"— gluon process. Note that 
the value of rric used is larger than that obtained in LO analyses of CC charm production. 
This choice tends to reduce the NC charm contribution to the WSM sample and results in 
more conservative limits on FCNC production. The NC charm quarks are fragmented and 
decayed using procedures adapted from the CC charm simulation, and the resulting WSM 
events are then simulated with the full MC. 



IV. RESULTS AND INTERPRETATION 



A. FCNC Production 



The neutrino FCNC u ^ c cross section can be parameterized to LO in QCD as 

da {Uf^u ^ Uf^c; c ^ H+) Kc ^ [ 2n, . 2 ni^ - v) - ^v/O] da {i^f,d ^ iJ,-c;c ^ jj,^ 

cos p + sm p- 



dC,dy 



cd 



1 - y + xy/i 



d^dy 

(4) 



Here Vcd is the c ^ d CKM matrix element, represents a possible u ^ c coupling, 
siv? [3 gives the fraction of right-handed coupling of the c— quark to the FCNC, y is the 
inelasticity, and ^ ~ a; (1 + ml/Q'^) is the fraction of the nucleon's momentum carried by 
the struck m— quark, with x the Bjorken scaling variable, the squared momentum trans- 
fer, and rric the effective charm quark mass. The d —>■ c charged current cross section 
da {v^d Ufj^c; c — > /d^dy is measured in the same experiment [?,?]. Since the u and 
d quark distributions are identical in an isoscalar target, the FCNC cross section should 
experience the same charm mass suppression as the analogous CC charm production. Frag- 
mentation of subsequent semi-muonic decays of charmed mesons should also be identical for 



We use the notation Vuc, Vbd, and Vsd in simple analogy to the CKM matrix in order to compare 
our results to those from FCNC decay searches.. We do not assume any constraints exist for this 
FCNC CKM-Uke matrix. In our notation, the FCNC left and right-handed couplings for charm 
are gl = |Kc|^ cos^ /? and gj^ = \ Vuc\^ sin^ p. 
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FCNC and CC-charm production. One therefore expects the extracted T^c to have httle 
model dependence. 

For FCNC bottom production there is as yet no measured CC analog final state. There- 
fore, the explicit LO QCD cross section, 

^ ' ^ = ^ K (1 - ^) (1 - ^^/^) + (1 - ^ + ^vio] (5) 

x(«(r,g^)+J(e',g^)), 

where M is the nucleon mass and E is the neutrino energy, must be convolved with h- 
quark fragmentation functions for mesons of type Bi {Di) and Bi meson decay distribution 
functions (A^) multiplied by appropriate branching fractions (F^) to yield a WSM cross 
section: 

da {ly^N vj)- b^fj,+) _ da (v^N u^bx) 

The struck quark momentum fraction ^' becomes ^' ~ x (1 + m^/Q'^) , with = 4.8 GeV/c^ 
the effective 6-quark mass. It is also possible for FCNC 6-production to form a WSM 
muon signal through the cascade 6 — > c — > This mode offers the advantages of the 
larger and higher ^ valence d-quark PDF at the cost of reduced acceptance for the softer 
c-decay muon. A similar expression holds for FCNC s ^ b transitions with the replacements 
u (^', Q^) + d (^', ^ 2s (^', g2), ^ \v,s\\ and sin^ /3' ^ sin^ 0' . 

Production cross sections for both c- and h- FCNC sources are computed from the 
GRV94LO PDF set [?] for several choices of right-left coupling admixtures. Acceptance 
for a charm FCNC- WSM signal is calculated using a fragmentation-decay model tuned to 
NuTeV and CCFR dimuon data [?]. For FCNC-WSM from 6-quarks, fragmentation and 
decays are handled with the Lund string fragmentation model [?]. Detector response is 
simulated with the full hit-level MC. 
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B. Fits to Data 



Binned likelihood fits are performed to the yvis distributions of the data using a model 
consisting of all conventional WSM sources described above and an FCNC source. The 
fit varies the level, but not the shape, of the FCNC signal contribution. The NC charm 
contribution is also varied in shape and level by allowing rric to fioat within its errors. The 
three FCNC sources {u ^ c, d ^ b, and s — > 6) are treated separately. Only neutrino data 
is used for the u ^ c, but both modes are used for FCNC bottom production to exploit 
the possibihty of a cascade decays to charm. A series of fits are performed for each FCNC 
source, corresponding to different mixtures of right and left-handed FCNC couplings to the 
quarks; a typical result is shown in Fig. ??. 

In all cases, the signal for FCNC is within ± 2.0 cr of zero, and hmits are set accordingly. 
Since Gaussian statistics apply, the 90% confidence level upper limit is set by adding 1.64(j 
to the best- fit value if the best-fit value is positive, or 1.64cr to zero if the best fit is negative. 
Here, a consists of the statistical error from the fit added in quadrature to the estimated 
systematic error described in the next section. Table ?? summarizes the fit results. 

C. Systematic Errors 

The dominant systematic errors result from modeling the rejection of CC charm events, 
and the overall normahzation of CC charm events. Estimates of systematic uncertainties 
are obtained by varying the event selection procedure as well as parameters characterizing 
the detector response and physics models. Errors are assumed to be independent. 

Charged current charm events are removed by requiring that exactly one track be found 
and reconstructed by the NuTeV tracking software. Another independent way to remove 
dimuons is to use calorimeter information. The stop parameter is the first of three con- 
secutive counters downstream of the interaction, each with less than 1.5 MIPs. The stop 
cut requires that the distance between the interaction and the stop counter be less than 15 
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counters. Replacing the tracking cut with the stop cut gives the systematic errors hsted in 
Table ??. 

The next largest systematic error is due to the normalization of CC charm events. Nor- 
malization of these events is obtained from the right-sign muon CC sample. One can also 
normalize CC charm events with only one reconstructed track, to those with both tracks 
found. These normalizations disagree by 3% resulting in the systematic errors listed in Ta- 
ble ??. Systematic errors due to the beam normalization, detector calibration, and other 
sources are small. 



D. Comparison to Limits from Decays 



For comparison purposes, the following expressions are used to relate FCNC heavy flavor 
meson decay branching fractions (BF) to the parameter T4c- 



BF (d° ^ £+r) = 2 

BF (d+ TT+tr) 
BF {p^ K+t£-) 



^ ml 



V ^ 



BF [d^s ^ /^""^J , 
BF (l>+ ^ TT+^+i/^) , 
BF (d^ ^ r^tv,) . 



For estimates of VcU) and Vgh from B decays, it is assumed that 

Vm 



2 2 

mi 



BF (B^ in-) = 2 ^ ^BF {B- ^ ,n,) , 

BF (b+ ^ 7r+£+r ) - ^ ' (5° ^ 7r-£+Ui) , 
BF (5° ^ in-) = 2 ^ ' ^BF (B^ ^ ^n,) , 

BF (b+ ^ K+e+r) = Y!^\f{b+ ^ DH+ut) . 



(7) 
(8) 
(9) 



(10) 

(11) 
(12) 

(13) 

Measured values [?] are used for the branching fractions on the right hand side except for 
the leptonic decay B^ —>■ /I'^Vf^, for which it is assumed that 



BF (b+ /x+i/^) = 2.2 X 10"^ (/b/200 MeV)^ , 



(14) 



11 



with fs — 200 MeV, the B decay constant. 

1 2 I 1 2 1 2 

Table ?? summarizes the hmits on \Vuc\ , , and \Vsb\ from meson decays. We note 
that our overall limits from neutrino scattering, which would approximately correspond to 
decay searches of the type D — > ly^Pf^X and B — > i/^P^X, are generally weaker than the 
decay search limits. Our result for Vdb is competitive, and we have effectively added new 
modes to the search that do not depend on specific mechanisms for heavy meson decay. 

V. CONCLUSION 

In this paper we have established a new method for probing FCNC processes in deep 
inelastic neutrino scattering. Our experiment tests for FCNC at the inclusive quark level, 
and we are particularly sensitive to any FCNC process in which the mediating field couples 
more strongly to neutrinos than to charged leptons. We observe no evidence for FCNC 
interactions, and we set limits on the effective mixing elements |Kic|^) and iHsl^ at 

the 10-^ level. 
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FIGURES 



NEUTRINO MODE 




ANTINEUTRIINO MODE 




50 100 150 200 250 300 



50 100 150 200 250 300 



FIG. 1. Comparison of v and v CC energy spectrum for data (pluses) to MC (histogram) using 
the GEANT-based flux. 
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FIG. 2. Comparison of data to MC of Ehad/ {Ehad + Efj_2) for dimuon events with two 
toroid-analyzed muons. 
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FIG. 3. Comparison of yvis distributions of data (pluses) to predictions of all Standard Model 
sources (solid) of WSM's and the FCNC signal(dashed). The plot on the left is neutrino mode, 
while that on the right is anti-neutrino mode. 
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TABLES 



Source z^-mode(%) ^'-mode(%) 

Beam Impurity 67 83 

Charged Current Charm 19 8 

Charge Misidentification 5 5 

Neutral Current Charm 5 2 

Neutral Current -k/K decay 2 1 

Charged Current -k/K decay 1 1 
TABLE I. Percentage of WSM's for each source in a given mode. 

z^-mode i^-mode 

Ez/ > 20 GeV Ei^ > 20 GeV 

scraping 53% 24% 

charm 10% 25% 

KO 12% 16% 

other prompt 9% 22% 

muon decay 11% 11% 

K^TT^ H 5% 2% 



TABLE II. The percentage of beam impurities due to a given source in each mode. 
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(-1.6±0.90±0.6)xl0-3 


1.8x10-3 




0.90 


(-1.4±0.79±0.6)xl0-^ 


1.6x10-3 




1.00 


(-1.3±0.68±0.6)xl0-3 


1.5x10-3 


s^b 


0.0 


(-17.3±17.3±3.51)xl0-3 


29x10-3 




0.10 


(-13.6±6.6±3.3)xl0-3 


12x10-3 




0.35 


(-3.6±1.9±2.7)xl0-=^ 


5.4x10-3 




0.65 


(-1.9±1.0±2.0)xl0-^ 


3.7x10-3 




0.90 


(-1.4±0.7±1.5)xl0-^ 


2.7x10-3 




1.00 


(-1.3±0.7±1.2)xl0-3 


2.3x10-3 



TABLE III. Results of the FCNC fits. 
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Transition Coupling Dimuon Rejection Dimuon Normalization Energy 



Beam 



Total 



« — > c 


L 


0.30x10" 


-3 


0.32 XlO" 


-3 


0.23x10" 


-3 


0.01 XlO" 


'S 


0.50 XlO" 


-3 


u ^ c 


R 


8.08x10- 


-3 


3.22 XlO" 


-3 


0.24x10" 


-3 


0.02 XlO" 


-3 


8.70 xlO- 


-3 


d^b 


L 


0.36x10- 


-3 


0.22 xlO" 


-3 


0.51x10- 


-3 


0.14 xlO- 


-3 


0.68 xlO- 


-3 


d-^b 


R 


0.25x10" 


-3 


0.10 XlO" 


-3 


0.04x10" 


-3 


0.55 XlO" 


-3 


0.61 XlO" 


-3 


s^b 


L 


2.37x10- 


-3 


0.21 XlO" 


-3 


0.91x10" 


-3 


2.42 XlO" 


-2 


3.51 XlO" 


-2 


s^b 


R 


1.08x10" 


-3 


0.06 XlO" 


-3 


0.20x10" 


-3 


0.54 X ID- 


-2 


1.22 xlO" 


-2 



TABLE IV. Tabic of systematic errors on FCNC results. "L" refers to pure left-handed coupling 
(sin^ j3 = 0), while "R" refers to pure right-handed coupling (sin^ (3 = 1). 





FCNC 


BF 




Allowed 


\v\' 




Limit with 






decay 


Limit 




Dccay 


limit 




BR (UTor 


reference 






1.7 xlO- 


-5 




2.3 xlO- 


-4 


2.7 xlO-^ 


[?] 


S± 




3.9 XlO" 


-3 




1.6 XlO" 


-3 


2.1 xlO-3 


[?] 






3.9 xlO" 


-5 




2.4 xlO" 


-5 


2.1 xlO-^ 


[?] 



TABLE V. Limits on FCNC couplings from meson decay searches, with |Fp= |KtcP,|Vdf)P, or 
iFsbP as appropriate. Equations ??-?? relate branching fraction(BF) limits to the \V\'^ limits in 
the table. 
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